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On I Abstract 

Reconnection processes of twin-chains are systematically studied in the quadratic twist 
^ ' map. By using the reversibility and symmetry of the mapping, the location of the indicator 

points is theoretically determined in the phase space. The indicator points enable us to 
' obtain useful information about the reconnection processes and the transition to global 

chaos. We succeed in deriving the general conditions for the reconnection thresholds. In 
addition, a new type of reconnection process which generates shearless curves is studied. 

> ' 

Csl ' 1 Introduction 
O 

i In the past decades, enormous effort has been dedicated to the study of two-dimensional area- 

00 ' preserving maps with the twist condition]^], but very few studies have been made of nontwist 

maps. Recent studies on nontwist maps have revealed that rich properties are generated by 
violating the twist condition. [2-6] 

In a previous paper we studied the properties of the quadratic twist map and numerically 
determined the critical boundary in the two-dimensional parameter space, where the transition 



' to global chaos occurs. The critical boundary has many sharp singular structures, and their 

1^ ■ locations seem to have a one-to-one correspondence with those of the reconnection thresholds. 

The relationship between the transition to global chaos and the reconnection processes was first 
pointed out by Howard and HohsQ, but it has not yet been thoroughly investigated. 
^ ' In order to investigate the detailed structure of the critical boundary, one needs accurate 

■ information regarding the reconnection processes. In this paper, we study the details of the 

reconnection processes in the quadratic twist map and propose a theoretical method to deter- 
mine the reconnection thresholds. We show that the reversibility and symmetry of the mapping 
guarantee the existence of the "indicator points" in the phase space. These enable us to study 
the reconnection processes systematically. For twin-chains of period one and period two, the 
reconnection thresholds have already been determined, either exactly or approximately |P, |3|. 
The method presented here reproduces results which have been previously obtained, and it 
provides general conditions for the reconnection thresholds. 
The quadratic twist map(QTM) is defined by 

• \ In+i = In-Ksm{9n), ^ ■ ' 

f^{I) = 27Tfl - I\ (1 • 2) 
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where K represents the strength of the perturbation and fi the maximum value of the twist 
function f^{I). As the variable 6 is 27r-periodic, /i is also periodic with period 1. Moreover, 
Eq. (1-1) is invariant under the transformation {K,I) i-^ {—K, —I). Thus it is sufficient for us 
to consider the parameter regions n € [—0.5,0.5) and K e [0, oo]. 

The mapping Eq. (IT) is often called the logistic twist map or the standard nontwist map, 
depending on the form of the twist function /^(I).[2-5] In this paper, we adopt a form of 
Eq. (1-2) for which the twist condition 

^^/O for^I (1-3) 

fails at / = for any value of fx in the integrable limit {K = 0) and for which the parameter fi 
itself represents the maximum value of ffj,{I)- 

This paper is organized as follows. In §2, we derive the location of the indicator points which 
play an important role in studying the reconnection processes in the QTM and give a review of 
the previous results on the transition to global chaos. In §3, the reconnection processes of even- 
periodic twin-chains are studied by using the indicator points. The reconnection thresholds are 
analytically derived, and a new type of reconnection process is analyzed. In §4, we focus on the 
reconnection processes of odd-periodic twin-chains, and a new numerical method to determine 
the reconnection thresholds is proposed. Section 5 contains a summary and discussion. 

2 The indicator points of the QTM 

By reversibility, the QTM can be rewritten as T = M2M1, where Mi and M2 are given by 

M.:[';,Z-'^''^-'''' (2.2) 

and they satisfy = Ml = 1- Moreover, the mapping T commutes with the mapping S (i.e., 
ST = TS) defined by 




(2-3) 

Previously!^, we showed that when there exists only one shearless curve in the phase space, it 
is an invariant set of both the mapping T and the mapping S. As shown in Appendix A, the 
invariant sets are given by the following sets Ij [j = 1,2): 

00 

x, = U {r"xf)}-o, (i = i,2) (2-4) 



m=— 00 



where Xi™^ and x"^'^ are the solutions of 

Mix = K^Sx (2 • 5) 

and 

M2X = K^Sx, (2 • 6) 
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respectively. Here m represents an arbitrary integer, and R represents the linear transformation 
defined by 



R 



The solutions of Eqs. (2-5) and (2-6) are as follows: 
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(2-9) 



Note that a;^™^ depends only on while x^^^ depends only on /x. If we use the 27r-periodicity 
of the variable Q in the interval (— 7r,7r], it is sufficient to consider only two terms (m = and 
m = 1) in Eq. (2-4). We shall call these points, x^™'' (j = 1,2), the "indicator points" in the 
present paper. As will be shown in what follows, we can detect the occurrence of the transition 
to global chaos and of the reconnection of twin-chains by investigating the iterates of the indi- 
cator points. 

When the iterates of the indicator points x^™^ are confined in a bounded region, the follow- 
ing two cases are possible: i) they are confined on a shearless curve (see Fig. 1(a)), or ii) they 
are not confined on a certain KAM curve, but several robust KAM curves surrounding them 
prevent global chaotic motion (see Fig. 1(b)). In either case, one can say that the bounded 
motion of the iterates of x^"^ ensures the existence of KAM curves. On the other hand, the 
unbounded motion of the iterates of x^^^ guarantees the non-existence of KAM curves. 

As reported in Rcf. 6), we have numerically determined whether the iterates of x^™^ are 
bounded or not for each set of ii and K. The phase diagram so obtained is shown in Fig. 2, 
where the iterates of x^"^ are bounded in the gray region, but unbounded in the white region. 
Thus, the boundary between the gray and the white regions displayed in Fig. 2 determines the 
critical boundary, where the transition to global chaos occurs. In numerical calculations, wc 
consider iterates to be bounded when the absolute values of their / components do not exceed 
2 during 10^ steps of iteration (i.e., < 2 for ^ra < 10^). These conditions are sufficient to 
detect the bounded motion, because the robust KAM curves are always localized around the 
indicator points x^"*\ as shown in Figs. 1(a) and (b). 



3 Annihilation and reconnection of even-periodic twin-chains 

In this section, we consider the case where either x^^^ or x^^^ is an even-periodic point with 
rotation number P/Q, i.e., 

rPt^xY' = 4™ ) , (j = 1, 2; m € Z) (3-1) 

where P and Q are relatively prime integers and let Q be even. By the symmetry of the 
mapping, Eq. (3-1) is reduced to more simple forms (see Appendix B): 

T^/^xS"^) = x^^^''^ (3 • 2) 



3 



and 

r-«/24-) = 4— (3.3) 

By solving Eqs. (3-2) and (3-3), one can obtain two independent relations between n and K 
parametrized by P and Q. Note that these relations do not depend on m at all. The exact 
iterates of x^^^ and X2^^ are given for Q = 2,4,6 and 8 in Appendix C. The obtained fx-K 
relations are summarized in Table I, where we put 

= l + 2cos(27r/u- ^/C^) (3.4) 

and 

, ^ cos(37ru - cos^fvr/i)) ,^ ^, 

r//i,K =1 + '-^ (3-5) 

cos(7r/ij 

Let us denote the fi-K relations derived from Eqs. (3-2) and (3-3) by e[^^^'^ and 
respectively; i.e., 

i?f/«) = {(/.,i^) |fp,Qxf)=xf)}, (j = l,2) (3-6) 
where m is an arbitrary integer and 

fp,Q = R^T^. (3 • 7) 

In the \i-K parameter space, the two sets e'^^^'^^ and e'^^^^ are described by two different 
curves which pass through the trivial point (|U, K) = {P/Q, 0). These curves are shown by solid 
lines in Fig. 2 for several values of P/Q (those for Q < 8). There are two types of characteristic 
fj,-K curves given for each value of P/Q- We shall call these curves the right {P/Q)-curve or the 
left {P/Q)-cuTve, depending on their relative positions, corresponding to the right-hand side or 
the left-hand side in Fig. 2. For (P/(5)-curves that have a crossing point, we only apply the 
notion of right and left below the crossing point. Some of the right and the left {P/Q)-cuvves 
are summarized in Table II. 

Note that the {P/Q)-curves are quite similar to the (P/(5)-bifurcation curve introduced 
by del-Castillo-Negrete et al.[Q], where the (P/Q)-bifurcation curve is defined as the fi-K locus 
which represents the creation/annihilation points of the twin-chains with rotation number P/Q. 
It will be shown later that either the right or the left (P/Q)-curve which we have introduced 
here coincides well with the P/Q-bifurcation curve. 

From numerical calculations, we can find the following general characteristics of the iterates 
of the indicator points x^"^^ and the characteristic Ej^^^^ curves defined by Eq. (3-6): 

(i) For the fi-K parameter values given by e[^^^\oi £^2^^*^^), successive iterates of X2^\or 

x^^) converge to a periodic point with rotation number P/Q under the mapping Tp^Q. 

(ii) For the ^-K parameter region between the right and the left {P/Q)-curves, successive 

iterates of both Xi^^ and X2^^ converge to an unstable periodic point with rotation number 
P/ Q under the mapping Tp^q . 

These properties are numerically illustrated in Fig. 3 for the case of P/Q = —1/2. 

It is found in many cases that twin-chains with rotation number P/Q (hereafter referred to 
as (P/Q)-twin-chains) are created or annihilated at parameter values on the left {P/Q)-cuive, 
and that the reconnection of the (P/(5)-twin-chains takes place at parameter values on the 
right (P/(5)-curve. Thus, the left {P/Q)-cmve and the right (P/(5)-curve correspond to the 
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annihilation threshold and the reconnection threshold, respectively. In Fig. 2, we can see nice 
agreement of the critical boundary and the right (P/(5)-curves for many P/Q. This is more 
clearly shown in Fig. 4 for the case of P/Q = —1/4. However, this correspondence does not 
always hold for all P/Q. For instance, in the cases of P/Q = —1/6 and P/Q = —1/8, it is the 
left {P/Q)-curve that coincides with the critical boundary. This is more clearly shown in Fig. 5 
for the case of P/Q = —1/8. This exceptional correspondence occurs when the left (P/Q)-curve 
changes to correspond to the reconnection threshold at a certain point on the curve. 

Let us illustrate this for the left ( — l/8)-curve. Figures 6(a)~(e) display phase space por- 
traits for various values of K, where K increases along the left (— l/8)-curve. At K = 0.5, 
marginally stable periodic points with rotation number —1/8 exist in the phase space, as shown 
in Fig.6(a). At K = Ki 0.519, the periodic point becomes unstable and gives birth to two 
elliptic points, as shown in Fig. 6(b). We refer to the chain of these periodic points as the 
primary chain. At K = K2 — 0.5492, another pair of periodic points with rotation number 
—1/8 is created on both sides of the primary chain via a saddle- node(SN) bifurcation. Figure 
6(c) exhibits the phase space portrait at if = 0.552, where two chains of SN pairs appear, in 
addition to the primary chain mentioned above. We refer to these chains of SN pairs as the 
secondary chains. Moreover, at K = ~ 0.561, the two separatrices of the primary chain 
and those of the secondary chains merge, as shown in Fig. 6(d). As a result of the reconnection 
process induced by the merging of the primary and secondary chains, vortex pairs and two 
Poincare-Birkhoff chains come to appear, as shown in Fig. 6(e). 

To summarize, the left (— l/8)-curve corresponds to the annihilation threshold of the pri- 
mary chain for K < Ki. While, for K > 1^3, it corresponds to the reconnection threshold of the 
vortex pairs that are formed via the reconnection process of the primary and secondary chains. 



4 Reconnection of odd-periodic twin-chains 

In this section we characterize the reconnection of odd-periodic twin-chains by using the indi- 
cator points. At a reconnection threshold of the (P/ (5)-twin-chains, successive iterates of both 
^(m+P) ^___^(m) g^ppj.Qg^g]j ^]jg game hyperbolic periodic point of the reconnecting twin- 
chains under the mapping Tp^g. Figure 7(a) corresponds to the case of the (— l/3)-twin-chains 
at the reconnection threshold. Figure 7(b) is the magnification of the phase space near the 
lower hyperbolic periodic point of the (— l/3)-twin-chains, where xf*^ and Tl'i 3 T^X2'^ 

approach the same hyperbolic periodic point as n increases. We have numerically confirmed 
this property for all odd-periodic twin-chains with period Q < 9. This property seems to hold 
even at a large ii'-value, as shown in Fig. 8, where we can still observe monotonic convergence 
in spite of strong chaos around a hyperbolic periodic point. 
Consider the quantity dXn given by 

sx^ = II - f^,Q r^x(-)||, (4 . 1) 

where m is an arbitrary integer and 1 1 • 1 1 represents the Euclidian norm. At the reconnection 
threshold, SXn is expected to decrease exponentially as 

SXn oc A", (4 • 2) 

where A represents the smaller (stable) eigenvalue of the tangent map of the mapping T'^ 
evaluated at the hyperbolic periodic point of the (P/(5)-twin-chains. Numerical evidence for the 



5 



exponential scaling regime of Eq. (4-2) at the reconnection threshold of the (— l/3)-twin-chains 
is shown in Fig. 9 for various K-values. There appear large deviations from the scaling regime 
when the number of iterations n becomes extremely large, because the monotonic convergence 
of 6Xn is violated by heteroclinic chaos near the hyperbolic periodic point. However, as shown 
in Table III, the decay rates in the scaling regime are in good agreement with the theoretically 
estimated values of ln(A). 

On the basis of the remarkable convergent property mentioned above, the determination 
of the reconnection threshold is reduced to the problem of finding the parameter values where 
< e holds for a given small value e with a certain large A^. The merit of this computational 
criterion is that it does not require any additional knowledge about the detailed structure of 
twin-chains, such as the location of the periodic point and its stability. Applying this criterion, 
we numerically determined the reconnection thresholds for odd-periodic twin-chains with the 
period Q <9, which are indicated by solid lines for Q=odd in Fig. 2. The results based on the 
above technique seem to reproduce the reconnection thresholds. This is illustrated more clearly 
in Fig. 10 for the case of the (— l/3)-twin-chains. 

5 Summary and discussion 

We have studied in detail the reconnection processes of twin-chains in the QTM by investigating 
the iterates of the indicator points Xj^^ which belong to the invariant sets of both the mapping 
T and the mapping S. For the case of even-periodic twin-chains, we introduced {P/Q)-curves 
in the parameter space, on which either the annihilation or the reconnection of the {P/Q)-tw'm- 
chains occurs. We succeeded in deriving the {P/Q)-cuTves analytically for Q < 8. One can 
also derive {P/Q)-curves with larger value of Q by solving Eqs. (3-2) and (3.3). On the other 
hand, for the case of odd-periodic twin-chains, we propose a numerical method to determine 
the reconnection threshold. The method is based on the observation that the successive iterates 
of the indicator points converge to the same hyperbolic periodic point at the reconnection 
threshold. Numerical results clearly show that the critical boundary is described well by the 
reconnection thresholds of twin-chains. 

Finally, we briefly discuss the reconnection process between the primary and the secondary 
twin-chains which was found for the (— l/8)-twin-chains. The variety of reconnecting twin- 
chains has been studied in cubic and quartic twist mapsj^. For example, cubic twist maps 
exhibit the reconnection of three island chains with the same rotation number. This is a 
natural consequence of the fact that the cubic twist function generically has two extrema and 
that each extremum generates an inherent reconnection process. However, it is surprising that 
multiple reconnection processes occur simultaneously even in the QTM whose twist function 
has only one extremum. The occurrence of the secondary twin-chains inevitably ensures the 
existence of multiple shearless curves in the phase space, neither of which passes through the 
indicator points Xj^\see Figs. 6(c) and (e)). Our result suggests that shearless curves can be 
created by nonintegrable nonlinear perturbations. This problem will be discussed elsewhere. 
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Appendix A The Invariant Sets X, (j = 1, 2) 

It is trivial to show that Ij (j = 1, 2) given by Eq. (2-4) is an invariant set of T. Here we show 
that Tj is invariant under the transformation 5, i.e, SXj = Xj (j = 1,2). In order to do this, we 
show that ST'^x'P^ is also an element of Ij as follows: 

= r"-iM2xf "'^ = r"x(^--'™) , (A • 1) 

= T^+^M^xf"^^ = n^^-^""^ , (A • 2) 

where we have used the relation ST = TS in the first equalities, and Equations (2-5) and (2-6) 
in the second equalities. Moreover, we use the following relations in the last equalities: 

M2xi™) = rxi'-"\ (A -3) 

M,x^^^ =T-'x^^-"'\ (A -4) 



Appendix B Derivation of Eqs.(3-2) and (3-3) 

Assume that Xj^^ {j = 1,2; 
rotation number P/Q; i.e.. 



Assume that Xj^^ {j = l,2;m G Z) given by Eqs. (2-8) and (2-9) are even-periodic points with 



The property of the involutions Mi and M2 implies 



i2^r«xf)=arfMi = l,2) (B • 1) 

Ml and M2 implies 
^ T-^^Mi, MsT-'^/^ = T'^/^Ms. (B • 2) 

Therefore we have 

MiT^'^^r^ = T-^/^Mixi'") = T-^'''W'Sx^r^ = R^'ST-'^/^x^;^^ (B • 3) 

and 

M2r-«/2arf ) = r«/2M2xf ) = T^/^R^Sx^^ = R^ST^'''xt^ , (B • 4) 

where we have used Eqs. (2-5) and (2-6) in the second equalities and the relation ST = TS in 
the third equalities. 

On the other hand, Eq. (B-1) implies 

r-Q/2^(-) ^ RPTQI^xf\ (j = 1, 2) (B • 5) 

rQ/2^M^^-Pr-Q/2^M_ (_^- = i,2) (B-6) 

Substituting (B-5) and (B-6) into (B-3) and (B-4), respectively, we have 

and 

MsT-'^/^a-M = Rr^-PsT-'^/^xt^ . (B • 8) 

Since these equations are equivalent to Eqs. (2-5) and (2-6), the solutions are given by 

T«/2xi™) = (B.9) 

r-Q/2^M^^(m-P)^ (B-IO) 

respectively. 
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Appendix C Exact Iterations of x^^^ and X2"^ 

The exact expressions for T'^x^^^ (n = 0, ±1, ±2) are given by 



T-'^x 



Tx 



T'^x 



rn) 
m) 
rn) 
m) 
m) 



{2m-l)-AT:^i + \K\l + i{^i,Kf) 



(-l)-+iii^{e(;u, K) + 2 cos(47rM - \K^{1 + ^{li, KY))] 
f(2m-l)-27r^ + ii^2- 

f(2m-l)' 

(_l)m+ll^^ 

f (2m- 1) + 27r/x- 

(-l)-iK 
f (2m - 1) + 47r/x - ^^^(l + ^(/x, if)^)' 
(-l)-iifC(M,if) 

where ^{iJ,,K) is given by Eq. (3-4). 

While the exact expressions for T^'x'^^ (n = 0, ±1, ±2) are given by 



T-'x. 



(m) 



f (2m - 1) - SiTfi + i^^ cos^(7r^) 
(-l)m+i K cos{7r n)r](ii, K) 

f(2m- 1) -vr^ 
(_l)m+i^cos(7r^) 

f (2m- 1) +7r//' 


f (2m - 1) + 2,71 II - K'^ cos^inn) 
(-l)™i^cos(7r^) 

f (2m - 1) + - cos2(7r//)(l + ry(/x, K)^) 
K cos{Tr iJ,)ri{tJ,, K) 



(C-1) 
(C-2) 
(C-3) 
(C-4) 
(C-5) 



(C-6) 
(C-7) 
(C-8) 
(C-9) 
(C-10) 



where r]{ii,K) is given by Eq.(3-5). 
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Tables 



Table I. The iJ,-K relations derived from Eqs. (3-2) and (3-3). 



Q 


P 


The 


K relation obtained from 


The IX 


-K relation obtained from 

RPTQxf^ = 


2 


-1 


8tt{h- 


- P/Q) = K' 


l^ = - 


1/2 


4 


±1 


Sirlfi - 


- P/Q) = 


47r{fi - 


-P/Q) =K'^ cos2(7r^) 


6 


±1 


247r(// 


-P/Q) = K\2 + ^{f,,K)^) 


37r(;U - 


-P/Q) = K'^ cos2(7r^) 


8 


±3,±1 


167r(/x 


-P/Q) = K\l + i{^i,Kf) 


87r(/n- 


- P/Q) = cos2(7rAi)(2 + 77(m, 



Table II. The relative location of the two curves 
corresponding to e[^^^^ and £2^^^^ . 



Q 


P 


left 


right 


2 


-1 


E2 


El 


4 


-1 


El 


E2 




+1 


E2 


El 


6 


-1 


E2 


El 




+1 


E2 


El 


8 


-3 


El 


E2 




-1 


El 


E2 




+1 


E2 


El 




+3 


E2 


El 



Table III. Slopes of the scaling regime in Fig.9 
and theoretically estimated values of ln(A). 



K 


slope 


ln(A) 


0.40 


-0.140 


-0.142 


0.60 


-0.320 


-0.325 


0.85 


-0.640 


-0.668 
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Figure captions 

Fig. 1. Two types of bounded motions, (a) {fi,K) = (-0.366328694,0.80), (b) = 
(0.089259, 0.28). In both figures, 50 successive iterates of x^^^ and x^2^ under the mapping 
T are plotted by sohd circles (•) and crosses (x), respectively. 

Fig. 2. Phase diagram of the quadratic twist map. The iterates of the indicator points x^™"^ are 
bounded in the gray region but unbounded in the white region. The boundary between 
the gray and the white regions corresponds to the critical boundary, where the transition 
to global chaos occurs. The solid lines represent the (P/Q)-cutvcs for cvcn-pcriodic twin- 
chains with period Q < 8 and the reconnection thresholds for odd-periodic twin-chains 
with period Q <9. 

Fig. 3. Phase space portraits and successive iterates of x^^ and x^^ under the mapping r_i^2 = 
R~^T^, where solid circles (•) represent the iterates oi Xi^\ and crosses (x) represent those 
of X2^\ (a) At {i-i, K) = (—0.5,0.5) G E2^^^\ x^2^ is a marginally stable periodic point, 
and successive iterates of x^^^ converge to x^"^ under the mapping r_i,2- (b) At {fi,K) = 
(-0.490052816,0.5) £ e[^^'^\ x^^'* is an unstable periodic point, and successive iterates of 
^2^^ converge to x^'^ under the mapping T_i,2- (c) At {lJi,K) = (—0.495,0.5), successive 
iterates of both x^'^ and X2^^ converge to an unstable periodic point under the mapping 

f_l,2. 

Fig. 4. Magnification of the phase diagram and the right and left (— l/4)-curves. 

Fig. 5. Magnification of the phase diagram and the right and left (— l/8)-curvcs. 

Fig. 6. The reconnection process of the primary and secondary chains with rotation number 
P/Q=— 1/8. The values of /j, and K increase along the left (— l/8)-curve as follows: (a) 
{fi,K) = (-0.084975,0.50), (b) {n,K) = (-0.079063,0.53), (c) (^,K) = (-0.074355,0.552), 
(d) {ii,K) = (-0.072334,0.561), (e) {ii,K) = (-0.070258,0.57). 

Fig. 7. (a) Phase space portraits of the (— l/3)-twin-chains at a reconnection threshold {fi,K) 
= (-0.32681945,0.4). (b) Magnification of the box in Fig. 7(a). 90 successive iterates of 
xf*^ and T X2'^ under the mapping ^-1,3 axe also shown as solid circles (•) and crosses (x), 
respectively. 

Fig. 8. (a) Phase space portraits of the (— l/3)-twin-chains at a reconnection threshold (jU,if) 
= (—0.30159859,0.85). (b) Magnification of the box in Fig. 8(a). 12 successive iterates of 
x^'^^ and T x^^ under the mapping T^i,3 are also shown as solid circles (•) and crosses (x), 
respectively. 

Fig. 9. Numerical evidence for the exponential decrease of SXn at a reconnection threshold for 
K = 0.4, K = 0.6 and K = 0.85. 

Fig. 10. Magnification of the phase diagram and the reconnection threshold for the (—1/3)- 
twin-chains. 
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